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Abstract. We show that the Friedmann-Lemaitre-Robertson- Walker equa- 
tions with scalar field and perfect fluid matter source are equivalent to a 
suitable non-linear Schrodinger type equation. This provides for an alter- 
nate method of obtaining exact solutions of the Einstein field equations 
for a homogeneous, isotropic universe. 



1. Introduction 

Recently, there have been interesting reformulations of Einstein field equa- 
tions for scalar field cosmologies (both for isotropic and anisotropic mod- 
els) in terms of generalized types of Ermakov-Milne-Pinney (EMP) equations; 
see PP, ini, |H]> for example. Such equations occur in a variety of physical 
contexts and in particular have served to provide a link between gravitational 
and non-gravitational systems |5] . We present in this paper an alternate (non- 
EMP) formulation of homogeneous, isotropic scalar field cosmology. Namely, 
we provide a formulation in terms of a non-linear Schrodinger type equation. 
Some applications to exact field solutions are presented, including some string- 
inspired cosmological solutions. 

One can set up a direct correspondence between EMP solutions and Schrodinger 
solutions which suggests, for example, that possible connections between this 
work and that in 5 can be pursued. The Schrodinger equation here also seems 
to be of some independent interest. 



2. Einstein Equations 

The Einstein equations for a Friedmann-Lemaitre-Robertson- Walker (FLRW) 
homogeneous, isotropic universe with scalar field 0, potential V, and perfect 
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fluid matter source assume the familiar form (for a vanishing cosmological con- 
stant) 

(2.1) H + — = -—PT , - = — 5— [PT + 3pt) , 

where the total energy density pT and pressure pt with matter contributions 
Pm, Pm are given by 



(2.2) 


Pt = P4, 


+ Pm, Pt 


for 






(2.3) 


P0 = Y + ^ ° <p 




(2.4) 






a" 




with D, n = 


= constants > 0, n 


' ' a 


scale factor 


a = a(i), and /c — 


0, 1, or -1 = 



2 

(n- S)!) 



the curvature parameter. G is 
Newton's constant and units are selected so that the speed of light is unity; 
cf. 0. Equations (i), (ii) imply the fluid conservation equation 

(2.5) pT + iH{pT+PT) = (), 

which by definitions (2.3), (2.4) reduces to the Klein-Gordon equation of motion 

(2.6) '<j) + ZH^ + V' o(j)^Q 

de f 

of the scalar field ip. Note that for 7j„ — ^ one has the equation of state 
= (7m - l)pm (by (2.4)). For 7„i = 1 (i.e. n=3), for example, Pm = 0, 
which is the case of a dust universe. 

There is also an equation of state p^ — (70 — l)p<;/, which follows by setting 

• r • 1 

70 — 2(jP' (jp' + 2[V o (f)) .70 however, unlike 7„i, is a non-constant function 
of time t. 

In this paper we set up a correspondence (a, 0, V) < — > u between solutions 
(a, (/), V) of the field equations (i), (ii) in (2.1) and solutions u of the following 
time-independent Schrodinger type equation 

(2.7) u"{x) + [E- P{x)] u{x) = -"^uix)^ 

with constant energy E and potential P{x). For a flat universe (fc = 0), or for 
the special values n=2 or 4 one notes that equation (2.7) is actually a linear 
Schrodinger equation. The correspondence provides an alternate method, for 
example, of solving the equations in (2.1). In some cases it is convenient to 
specify the scale factor a[t) a priori and then find the scalar field (p and potential 
V such that (a,(/), y) is a solution in (2.1). This can be done, in particular. 
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in our approach as a{t) is sufficient to determine u{x) in (2.7), which by the 
correspondence determines the desired (p and V. 



3. Description of The Correspondence (a, (p, V) < — > u 

With the above notation in place, and where we set = SttG for conve- 
nience, we can state the main theorem: 

Theorem 1. Let u{x) be a solution of equation (2.7), given E, P{x). Then 
a solution (a, (j), V) of the Einstein equations (i), (ii) in (2.1 ) can be constructed 
as follows, wliere D in (2.4) (which specifies Pt,Pt) is chosen to be "a : First 
choose functions cF(t), tp{x) such that 



(3.1) a(i) = u(a(i)) , i,'(xf 

Then one can take 



P(x). 



(3.2) 
(3.3) 



a{t) 
V 



u{<j{t))- 
12 



iu') 



/\2 



2u^P 12u^E Sku^ 



+ 



o V"^ 



Here, in fact, (a, (j), V) will also satisfy the equations 



(3.4) 



(3.5) V(m = ^ 



m 

H(t) 



a(tY 
2 , H(t) 



nP 
3a(t)" 

2k 



Za{tf 



+ 



in-6)D 
6a(t)" ' 



again D ^l]^ . Conversely, let (a,(j),V) be a solution of equations (i), (ii) 
in (2.1), with some D given in (2.4) that specifies pT, Pt in (2-2). Similar to 
(3.1), choose some solution a{t) of the equation 

(3.6) &{t)=a{t)-^. 

Then equation (2. 7) is satisfied for 

def. K^n^ 



(3.7) 

(3.8) 
(3.9) 



E 



P{x) 



12 



-D, 



u{x) 



def. 



a{a-\x))-^. 



Theorem 1 therefore provides for a concrete correspondence (a, (j),V) ^ u be- 
tween solutions {a,(p,V) of the gravitational field equations (i), (ii) in (2.1) 
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and solutions u ol the non-linear Schrodinger type equation (2.7). The so- 
lutions (a, (j), V) also correspond to solutions Y of the generalized Ermakov- 
Milne-Pinney equation 

(3.10) Y" + QY = ^ 



y3 ' 

as discussed in pi; also see [J, In turn, one can set up a correspondence 
Y ^ u between solutions Y of (3.10) and solutions u of (2.7), and thus obtain 
the correspondence {a,(j),V) <-> u - and the proof of Theorem 1. However, 
we prefer to give a direct proof of Theorem 1 - one that does not rely on the 
results in nor on the correspondence Y ^ u - although the latter route 
served as the motivation for the formulation of equation (2.7), and consequently 
of Theorem 1. It is implicitly assumed in definition (3.3) that the inverse 
function of -0 exists, which would not be the case if P{x) — 0. Therefore 
we first assume that P{x) is not the zero function. The case P{x) = will be 
discussed later. 

For the proof of Theorem 1 we first establish the salient formulas (3.4), (3.5), 

def n 

given E, P{x), u{x) in (2.7). By (3.2), uo a ~' which differentiated in 
conjunction with (3.1) gives (u'ocr)(uo(T) = —^a^^~^a —^{uoa)H. That 
is, u' o a = —jH which we also differentiate to obtain {u" o a){uo a) — —^H 
(again by (3.1)); i.e., 

n Ti Ti • 

(3.11) uoa = a~^, u'oa — ——H, [u" o a){u o a) = ——H . 
From equation (2.7), Pu^ = u"u + Eu^ -f so that by (3.11) 

(3.12) (Pu2)oa = -^ij-fi;a-" + ^a-2. 
Using (3.11), (3.12) we now obtain by definitions (3.2), (3.3), 

T. ^ 12 / 71 \2 2 f n- ^ _„ nk 12E _„ 3fc _2 
Voch = ^ ^ H] I --H + Ea "H a ^ H "H ^ 

2,H^ H D (11 \ 2k 2 , • def. 12E ^ 
= -F7T + ^ - 1 + 77?a smceZ) = ^— r), 

which is equation (3.5) as desired. Again by (3.1), (3.2), (p = o a){uo(j) ^ 



-^H - if-a-'' + j^a-^, which is equation (3.4). 



A NON-LINEAR SCHRODINGER TYPE FORMULATION OF FLRW SCALAR FIELD COSMOLOGY 



With equations (3.4), (3.5) now established, we can compute pr, Pt in (2.2), 
using (2.3): 



+ 



H 

(n-6) 



o 



3 ' 



Pt = 



dg. 
dg. 
(i)' 



g -£>a-" (by(2.3), (3.4), (3.5)) 
(3fca-^ + 3H'^)K-'^ - Da-"- 
{3ka-^ + 3H^)K-^ - pm (by(2.4)) 



P4> + Pm 

which is exactly equation (i) in (2.1). Similarly, by (2.3), (3.4), (3.5), 



P4> 



-{H-ka-^)K- 



D 6 6 D 



= -2HK-^ - ka-^R-^ - SH^K'^ + 



(3-n) 



Da-"" 



Pt 



dg. 

def. 



P<t> 



{2H + ka-'^ + 3H'^)K-'^ 

Prn = -i2H + ka-'^ + 



-p^ (by(2.4)) 



which with {i}' gives pT+SpT = 3{ka-'^+H^)K-^-3{2H+ka-^+3H^)K-^ = 
-6{H'^+H)K-'^ = -6|^^'-^ which is exactly equation (ii) in (2.1). The arrow 
u — > (a, (j), V) has therefore been established in one direction. 

For the other direction (a, (j), V) — > u, with (a, 0, V) a solution of the equa- 
tions in (2.1), we must show that u{x) defined in (3.9) solves equation (2.7), 
for the data E, P{x) defined in (3.7), (3.8), with cr(t) defined by (3.6). For 
convenience let g{x) = a~^{x) denote the inverse function of ait): a{g{x)) = 
X => &{g{x))g'{x) = 1. That is, 

(3.13) uix)g'ix) = 1 

since by (3.6), (3.9) this product is a{g{x))^^ g'{x) — (T{g{x))g' {x). Then 
u\x) = 



-^'^igix)) = ^a{g{x))g'{x) 



def. 



n 



'=^- --u{x)H{g{x))g'{x) 

71 

-^H{g{x)) 



u"{x) = --H{g{x))g'{x) 



n 
2 



" a(g(a^)) 
_a{g{x)) 



H{g{x)f 
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where by (2.1) the bracket here is 



K"^ k 
— g- [Pt{9{x)) + 3pt{9{x))] + ^(^g(^^^^2 - -3-pT(ff(a;)) 

4 2 

= — ^ [PT(ff(a;)) +pr(,9(.T))] +/cw(x)" (i.e. u{x)— ^' a{g{x))) 



2 r 



K 
2 



4>{9{x)f + 



nD 



M9{x)Y 



APlx) , 'nD , ,2 



+ ku{x)^ (by(2.2), (2.3), (2.4)) 

ku{x)^, by(3.8), (3.9). 



That is, we sec that u"{x) = -%9'{x) 



^^u{x) + ^u{x) u{x) + ku{x) 



P{x)u{x) + '-^Du{x) - ffcw(.T)^-i (by (3.13)), which by (3.7) is exactly 
equation (2.7). The proof of Theorem 1 is therefore complete, where we have 
assumed the existence of '^"^ in (3.3). 
Remarks. 



(1) The case P{x) = 0: 

In equation (2.7) we generally take a non-zero potential P{x). If P{x) = 
then il){x) is a constant function by (3.1) and therefore its inverse 
function 'tp~^{x) in definition (3.3) does not exist, which means that 
the expression for V there has no meaning. However, if P{x) = wc 
can instead define V{x) to be a constant function and we also define (/>(t) 
to be a constant function, because of (3.2). First note that equation 
(3.4) still holds (with the left hand side there being of course) since 



by (3.12), = -f ij - ^^Da- 
which multiplied by gives 



agam as £> = 



(3.14) 



2 • nD _„ 2k 



as claimed. Next note that the right hand side of (3.5) indeed is a 
constant function of t (for {P{x) = 0). Namely, differentiate (3.14) to 
obtain 



(3.15) 



H 



2k 
-— + 



6a" 



H 
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and then compute that 



d_ I _3_ 

di I 



H 



2k 



3 



6a" J 
(n - Q)nDH 



2HH 



-2k K^n'^D\ 4kH 
H 



3a2 



18a« 



= 3 



= 



2H 



2k 



3a2 

n^DH nDH 



(n - 6)nDH 



(by (3.15)) 



-nD 



6a" 

nDH 
a" 



6a" 
(by (3.14)) 



(as claimed) that for some constant Vo one has 



(3-16) 



2 , m) 



2k 



(n - 6)D 
6a(t)" 



= Vo. 



In summary, in case P{x) = in (2.7) we define (j){t) = any constant 
and V{x) = Vq in (3.16) (as definition (3.3) no longer has a meaning). 
Then equations (3.14), (3.16) (the proper versions of equations (3.4), 
(3.5)) hold. 

It follows that, again with a{t) u{a{t))^^ , where a{t) = u{a{t)), 
one does arrive at a solution {a,(j),V) of (i), (ii), in (2.1). Namely, in 
(2.2), (2.3), (2.4), p^{t) = Vo, p^{t) = -Vq, 



Pt 



= Vo + Da- 
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which is (i), and 



Pt 




— n 



=^ Pt + 3pT 



2Vq + (n - 2)Da 



—n 




{n - 6)D 



+ {n - 2)Da 



— n 



(by(3.16)) 




Ak 



+ 



2Dn 



3a" 




+ 




+ H 



6 a 



which is (ii). 
(2) Equations (3.4), (3.5) imply (i), (ii): 

The argument following the establishment of equations (3.4), (3.5) ac- 
tually shows (independently of u{x)) that if a{t) is given a priori, and 
if (t){t), V{x) are functions satisfying equations (3.4), (3.5), then auto- 
matically {a,(j),V) solves equations (i) (ii) in (2.1). 



As a simple illustration of the application of Theorem 1, choose u{x) — 
(1 - A'^uj'^x'^) jJ^ for A, w > 0, which solves equation (2.7) for E = ^, P{x) = 
2^2 (1 - A^u?) (1 - A^u?x'y^ = K^B^A^ (l - A^u;^x^y\ = 2{l~A-'uj^)/K^A^, 
n — 4, k ^ 1. In (3.1) we can take a{t) — ta.Bh{ujt)/A'^uj, ip{x) = ipo ± 

arcsin(A2ajx) — ± ^ arctan( "^Ji^-i^-i ) A^uj\x\ < 1 and obtain from 

(3.2), (3.3), a{t) — A cosh.{ujt) , (j){t) = V'oi^ arcsin(tanh(a;t)) — V'oi^ arctan(sinh(cji)) — 

■00 ± 2^ arctan(e'^*), V{x) = ^ + B"^ cos^ {^{x - '0o)), for constants ipQ, Vo> 
which is the EUis-Madsen solution in section 4.3 of 2 . We note that the 2 in 
the expression sin (2^((/) — 0o)) in equation (42) of (2j should not appear there. 
One can obtain in fact all of the solutions in P] for a suitable choice of u{x) 
and n. 

Ozer and Taha have considered in |S] two string- motivated solutions {aj ,4>j,Vj), 
j — 1,2, for fc = 1, Dj = 0, where the potentials Vj were specified a priori. 
One can also take the point of view of specifying the scale factors ai(t) — 
(oq -|- t^)2, a2{t) = flo -I- /2ao, for ao ^ 0, and then applying Remark 2 



to find {(jij.Vj). For a2{t), for example, H2{t) = 2t/(2a2 -t- t^), H^it) = 
2{2al - t^)/{2al + t^f , and one obtains by (3.2) (i>2{tf = At^/K'^{2al + 



4. Some Examples 
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2^ 



where we choose the positive square 



root. Also, H2it)^ + H2{tf/i + 2/802(^)2 = (lOt^ + 12ag)/8(2ag + t^f ^ 
V2{4>2{t)) = (lOt^ + \2al)/K'^{2al +t^f, by (3.5). For x > 4>o, we can take 

(t)2Hx) = V2\ao\ [e^i^'-M - l] ^ and compute that ^2(2:) = V2iM(l^2^ i^))) = 
|-5g-K(x-0o) _ 2e-2Jf(^-<^o)] 7X2^2^ Similarly, for the above scale factor ai{t) 



one can obtain via (3.4), (3.5) <j)i{t) = 0[,+-^log (l + ^yv{x) 
say for x > (f>Q. These solutions {aj,(f>j,Vj) can also be obtained by taking 
Uj{x) — aj{a^^{x))~^ (which is motivated by (3.2)) for the convenient choices 
rii = 4, ri2 = 2, and using the corresponding solutions (Ti (t) = a^^ a.Tcta.n{a^^t), 
-01 (x) = ipo - ;^ log(cos(aox)), (T2(i) = y2arctan((^/2ao)"^^), tp2{x) = ipo - 
log (cos (-^)) of the equations in (3.1) for Pi{x) = 4aQ tan2(aox), P2{x) — 
tan2(.^); El = E2 = 0. One can go beyond the assumption Di = D2 = 
(in |S]) and obtain solutions VjO(j)j via (3.4), (3.5). For example, for Di ^ 
one can check, using Mathematica for example, that 
(4.1) 



(jll{t)=<f) 

and that 



-arctanh 



DKH 



ao 



\/3iog {it + y^~3DK^ + m- 



(4.2 



ViiMt)) = 



3(4t2 



3a2) 



As a fourth example, and final one, for a parameter A G 
-ITStanh (a/^Ax 



.\ {0}, let u{x) = 

Then u(a;) solves equation (2.7) for _E = A + ^A2, 

v4 > 0, P{x) = A, n = 1, k = — 8A2. By an application of Theorem 1, one 
obtains the solutions (a, (f)±, V) given by 

16 



(4.3) 



a(<) - 
0±(O = 
y(a:) = 



1 + 6-1^/2"' ^ 

27 

± - A /^-^^ arcsinh 
3 V 3 Aif 



e 8%^^ 



37A2 



25/^2 

for integration constants c, 0o- 



135^ , 2 
^tanh 



3 /aAif, 

iV2Vl^^- 



00 ) 
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